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ABSTRACT ' 

The application of geometrical schemes similar to those of Buckminster 
Fuller's geodesic dome to large sphere antenna reflectors has been investi- 
gated, The purpose of these studies is to determine the shape and size of 
flat segmented surfaces which approximate general shells of revolution and 
in particular spherical and paraboloidal reflective surfaces. The extensive 
mathematical and computational geometry analyses of the reflector have re- 
suited in the development of a general purpose computer program. This pro- 
gram is capable of generating the complete design parameters of the dish and 
can meet stringent accuracy requirements. The computer program also includes 
a graphical "self contained" subroutine which enables one to graphically 
display the required design. 
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I. INTRODUCTION 


The last decade has witnessed a dramatic increase in the research activities 
dealing with the possibility of utilizing space for various commercial and 
scientific needs. In several recent issues of Astronautics and Aeronautics 
(see, for example, ,[1-4]), many articles have appeared which deal with diverse 
aspects of large space structures. These articles have identified various 
applications and also proposed novel designs of structures to meet such 
applications. A review of the research activities on space structures prior 
to 1966 has been documented in volume [5] that resulted from an International 
Conference on space structures. 

It has thus become necessary to find and analyze small lightweight 
structures that will be used easily to construct much larger space structures. 

It would be desirable for these structures to be isotropic in nature. However, 
construction requirements may make this infeasible, therefore requiring ortho- 
tropic or possibly completely anisotropic structures. The latter is also unde- 
sirable because of the added complexity to the problaa. Truss-type periodic 
(repetitive) structures have recently been analyzed' as candidates for 
space structures [5-8] . Here simplicity in construction coupled with large 
stiffness to density ratios will be most desirable. However, up to now, most 
of the extensive modeling and design of rod-like space structures have been 
concerned with flat structures in the form of either plates or three-dimensional 
Cartesian structures [7-9].-. 

Rod— like structures in- the form of spherical domes have previously been 
modeled and analyzed as candidates for many on-ground structural applications. 
These are knotm as the geodesic domes and constitute variations on the original 
spherical dome invented by Buckminster Fuller [10,11] . 
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The flat rod-like models studied so far in [7-9] are restricted in 
their possible applications to such structures as floors, ceilings and 
straight walls. However, the results reported in [7-9], although presented 
an initial stage in the understanding of the behavior of large space structures, 
gave insight and confidence into the possible extension to the study- of more 
geometrically complicated structures. Shallow shell structures have frequently 
been mentioned, for example, as candidates for building components in commun- 
ication systems, orbiting antenna, and solar energy satellites. 

In this report we describe methods of designing and estimating the geo- 
metric accuracy of rod-like shells of revolutions. These will consist of single 
surface shells in the form of either spherical or paraboloidal shells. Although 
our analysis will be kept general we do this with an important application in 
mind; namely the design and geometric analysis of large space antennas. 

Our ambitious aim of modeling general rod-like shallow shells of revolutions 
cannot be achieved in one shot and thus has to wait until we develop further 
basic understanding of some simplified special cases. Here we conceive of the 
general shell to also include the cases of concentric shells with rod-like in- 
ternal and external surfaces and cores., Having analyzed the flat rod-like 
structures the immediate difficulty in the study of the corresponding curved 
structures lies 'in the modeling Cspecification) of their geometry. In the 
case of the fiat structures the geometric arrangement is practically trivial 
since most of them can be constructed by using only a few ntimber of different 
length columns. As an example, we recall that the tetrahedral model of [7] 
and the- octetruss model of [9] use single rod length elements. 
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II. LITERATURE REVIEW OF SPHERICAL GEODESIC DOMES 


General Description ; 

The transition from a flat rod-like surface to its corresponding curved 
surface introduces troaendous amounts of geometric difficulties and constraints. 
Consider for example, the case of a (0, + 60°) flat plate arrangement shown in 
figure 1. Here ail rods have the same length. Now, in trying to deform figure 
1 to fit on a spherical surface for example, one finds it impossible to do so, 
especially if he insists on maintaining the lengths of the members to stay the 
same. In so trying one soon realizes that he does not have enough mathematical 
tools to construct such a discrete spherical shell. The geometric constraints 
encountered in constructing rod-like shells can only be understood and utilized 
if one familiarizes himself with the appropriate mathematical tools such as 
geodesics. Here geodesics is defined as the technique for constructing shell- 
like structures that hold themselves up without supporting columns . They could 
be very light and very strong. They can also be very large and hence attractive 
for space as well as earth structures. The geodesic spherical bubble erected 
to house the United States exhibits at EXPO |67 in Montreal is an example which 
exhibits all of the above attractive properties . Smaller geodesic domes have 
also been used as cabins, offices, playgrounds, and pavillians, etc. Yet, 
considering their apparent potential, in the quarter-century since Buckminster 
Fuller introduced them they have not been used very widely. This is perhaps due 
to the fact that they are mathematically derived structures and their mathematics 
has not been easily available. Parts for the self-supporting frame must be 
fabricated to close specification. The fabrication, with today's technology, 
is no problem, the problem is learning what the specifications should be. 
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FIG. 1. (0®, + 60°) ARRAY 
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In developing the geodesic dome, Tuller recognized that the equilateral 
triangle was the most basic geometric structure which is also inherently very 
stable and strong. By translating, rotating and piecing these triangles to- 
gether one can form regular polyhedra. One specif ic polyhedron, the icosahedron, 
was found to best approximate the sphere of all other polyhedronal forms . 

To see how this works, consider the twenty identical equilateral triangles 
of figure 2a. By cutting along the outer sides of this figure, wrapping it 
around and connecting the edges one gets the icosahedron of figure 2b. This 
icosahedron will have twenty identical triangular faces, twelve vertices and 
thirty equal lengthed edges. Moreover, there exists a unique sphere which 
circumscribes this icosahedron and passes through its twelve vertices. If one 
imagines blowing up the icosahedron to completely fill the sphere one recognizes 
that the twelve vertices maintain their original positions but the triangular 
sides become spherical and lie on the sphere. From the above discussions we 
can conclude that the icosahedron constitutes a rather rough approximation of 
a sphere. Better approximations can be obtained, however, by subdividing the 
individual icosahedron triangles into smaller ones and blowing them up to 
locate their vertices on the sphere. These subdivisions are known as the 
frequencies of the structure. The higher the frequency, the closer one gets 
to the sphere. Upon further subdivisions, one soon realizes, however, that the 
geometric constraints become enormous (as will be shown later, for example, the 
number of different lengths is' roughly in the order of the square of the 
frequency) . Various methods of subdividing the icosa triangles have been repoi^ted 
in the literature (see, for examples [10,11]). 

Simple geometric consideration will reveal that the side length L of the 
icosahedron is equal to 1.051462R where R is the radius of the sphere circum- 
scribing it. For further- discussion of the spherical geodesic dome descriptions 
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FIG. 2a 


ICOSA TRIANGLES 



FIG. 2b. ICOSAHEDRON 
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we'rrefer the reader to [10,11] . In what follows however, we shall treat 
the geodesic spherical dome as a very special case of our intended general 
modeling of rod-like spherical and paraboloidal caps . 

Having identified large space shallow antennas as important application 
of our rod-like shells we here point out the insufficiency and the inflex- 
ibility of the Fuller type geodesic spherical dome to meet their general 
geometric requirement. . For one thing in trying to extract a complete cap of 
the Fuller sphere it will be found out that such a cap ^■tIII have a unique focus 
to aperture diameter ratio (referred to as FOB) . Thus arbitrary degrees of 
cap shallowness will not be possible to model. Secondly, the number of rods 
emanating from the various vertices of the subdivision triangles (once on the 
sphere) will not be uniform. Specifically from each of the original ti^elve 
icosahedron vertices there ' emanate five rods where as from' all the remaining sub- 
division vertices there emanate six rod's. This nonuniformity will influence the 
^.degree of smoothness required. Thirdly, in constructing the Fuller sphere 
only the center of the sphere is used as the projection center of the sub- 
divisions. Fourthly, further development and analysis will be required to 
construct a Fuller type paraboloidal shell. The above unattractive properties 
of the Fuller shell will be more obvious later on in our modeling analysis once 
we introduce the many, nonrestrictive 'geometric degrees of freedom. 
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III. ANTENNAS IN THE EOEM OE ROD-LIKE SHELLS OE REVOLUTION 


In order to use rod-like structures as good candidates for building 
large space antennas we must first understand the global geometric description 
of the required structure. Specifically we must know if the antenna will be a 
part of a sphere, parabolid, ellipsoid or any other form of shells of revolutions. 

JU 

Once this is specified ^ two extra parameters such as the height, H and 

the aperture radius, R^, of the shell cup will he enough to completely specify 

the required geometry. 

A typical shallow shell of revolution cap is shown in figure 3a . By chang- 
ing the ratio of H to various degrees of shallowness will be realized. As 
for the rod-like approximation of the cap we proceed as follows : We subdivide 

the circumference of the aperture circle into n identical segments where n is 
an arbitrary integer. The points connecting these segments are numbered 2 - 
(n+1) and are then connected with thevertix point 1 (as shown for example in 
figure 3b for n=8) to form an n identical sided pyramid. In terms of H aud R^, the 
side lengths of the individual triangular face of 'the pyramid are given by 

1/2 

= (rJ + H^) (1) 

L 2 = 2R^ sin ^ (2) 

This pyramid will then form the roughest discrete approximation of the cap. 

Better approximations can be obtained, however, by subdividing the individual 
original pyramid triangular faces into smaller triangles and blowing them up 
to locate their vertices on the desired solid surface (see figures 4-6 for 
procedure illustration) . These subdivisions are 


If R^ and FOB are the given two parameters then H can be calculated as 


- {1 - 


H = 2 ROD X R. [1 - ^ 

^ leCEony 

a paraboloid. 


1/2 

-} ] for a sphere and H = R^/ (16 ROD) for 
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FIG. 4. PROCEDURE ILLUSTRATION OF SUBDIVISION AND 

BLOWING (n=6, N=7, FOD = 0.25) 
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PYRAMID 


PYRAMID 
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FINAL PARABOLOIDAL 
SHAPE 




FIG. 5. PROCEDURE ILLUSTRATION OF SUBDIVISION AND BLOWING 
(n=6, N=7, FOD = 0.25) 
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PYRAMID 



FINAL PARABOLOIDAL SHAPE 

FIG. 6. PROCEDURE ILLUSTRATION OF SUBDIVISION AND BLOWING 
(n=8, N=7, FOD = 0.25) 
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knovm as the frequencies of the structure. The higher the frequency is, the 
closer one gets to the required surface. Upon further subdivisions one soon 
realizes, however, that the geometric constaints become enormous. Such 
complexity of these constaints will be discussed later on in the analysis. 

III. 1 Pyramid* Face Breakdown' 

As we mentioned earlier the pyramid represents a rather rough approxi- 
mation of the required surface and that better approximations can be obtained 
by subdividing the individual face triangles into smaller ones and blowing them 
up so that their vertices will lie on the circumscribing surface. 

There are many types of breakdown for the original face triangle. We here 
mention two of them. The first is called the "alternate breakdown" in which one 
draws lines parallel to the sides of the triangle. The second hreakdotim is 
known as the "triacon" and it is obtained by drawing lines perpendicular to the 
triangle’s sides. Both breakdowns are illustrated in figure 7 for frequency 2. 




PIG. 7. TRIANGULAR BREAKDOWN (N=2) 
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There are many differences betoeen these two kinds of breakdowns; the most 
obvious two are; in the alternate breakdown the original triangle edges re- 
main part of the structure where, as in the triacon they do not and the second 
is that the alternate breakdown is possible in all frequencies whereas only 
even frequencies are possible in the case of the triacon breakdown. 

In the remaining of this report we shall concentrate on describing the 
alternate breakdown; the triacon breakdown will be discussed in a later report. 
Before we further discuss the form of subdivisions we shall first adopt the 
following appropriate coordinate system. 

The n-sided pyramid is oriented in the three dimensional rectangular co- 
ordinate system as shown in figure 8 . The origin of this system is dhosen to 
be the center of the sphere that circumscribes the p 3 rramid, namely the sphere 
that passes through the points 1-n. In terms of H and R^, the radius R of 
this unique sphere is given by 

R = (rJ + H^)/2H . (3) 

Accordingly, we choose the Z axis to pass through the vertex 1 with the X-Y 
coordinates being parallel to the aperture circle. Due to the symmetry of the 
pyramid faces we shall only treat a single one, namely the triangle 1,2,3 of 
figure 8. The projection of this triangle (1,2,3) on the horizontal plane 
circle is given by 1,2,3 and is shown together with the X-Y coordinate system 
in figure 9. With this choice of coordinate system the coordinates of the 
points 1,2,3 are given by 
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= (0,0,R) 

^ sin(^), R-H) 

= (\ cos(J), R^ slnCj), R-H) 

Notice that the coordinates of the point L’ are (0,0, R-H) . 


Accordingly, by subdividing the triangle 1,2,3 into frequency N as shown 
below in figure 10 we can determine the coordinates of the sub triangle vertices 

^IJ alternate subdivision rule (this rule has also been employed 

hy Glintoh flO] and others) . 


(X-X.) (X.-X-) 


(Y-Y^) CTo-Yj 

Y T T* I ■■ ^ + ‘ T ^ ^ 

JJ 1 ^ N N 


(Z -Z ) . (Z.-Z-) 

IJ ^1 ^ ^ N ^ ^ N 

where I and J are integers such Chat 

0 < I < I < N 


(5) 


ffatice from figure 10 that the distance between each two neighboring vertices 

L, 

along the 1,2' and 1,3 sides is constant and is equal to ^ and that the distance 
between each two neighboring vertices along the 2-3 side is a constant equal 
to , Now, since the coordinates of each vert^ are known, its distance D__ 
from the origin (0,0,0), -for example, is given- by 




1/2 


( 6 ) 


17 



1 



FIG. 10. BREAKDOWN NUMBEPaNG 
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III. 2 Projections on Spherical 'and Pa^boloidal ‘STtiTf aces 

In what follows we shall concentrate our efforts on developing discrete 
surfaces of spherical and paraboloidal cap surfaces. Simple geometric con- 
sideration will reveal that the equation of the sphere that passes through 
the vertices 1-n of the pyramid is given by. 


X + y + Z = R (7) 

As for theparaboloid that passes through the points 1-n, the appropriate 
equation is 


r3 

2^2 ^ 

X + Y = H" ^^"2) 


( 8 ) 


Csee sketch of figure 11 as an illustration of the paraboloid surface) , A 
schematic comparison of both the spherical and paraboloidal surfaces is shown 
in figure 12, 

III. 3 Pro.jection Centers 

- If the projection of the point (X^j, required surface 

is designated as ’ then the location of the points on the surface 

will depend upon their origin of projection. In what follows we shall leave 

the location of such a center arbitrary, namely (X^, Y^, Zp) . 

s s s 

The points (X^j, Y^^, can be obtained by connecting the points 

(Xp, Yp, Zp) and (X^j, Y^^, Z^j) by a straight line and extending it to 
intersect the required solid surface. Accordingly, the equation of this 
straight line is 


^ s „ V® V 7® 7 ^ a sphere 

^ij~^ _ ^ij~^p ^ ^ij~^p ® 

Xpj-Xp Xpj-Yp 2^j-Zp for- a paraboloid 


19 







FIG. 12. COMPARISON OF PARABOLIDAL AND SPHERICAL SURFACES 
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Once the specific surface (sphere or paraboloid) is specified - equation (9) 

s s s 

can be solved for individually in terms of the remaining 

quantities. For example, we obtain 

s s 

for a paraboloidal surface. - Similar expressions for and Z^j can be as 
easily obtained. Finally, substituting the resulting Y® and into 

the appropriate equation of required solid (equation (7) for the sphere and (8) 
for the paraboloid) one obtains the following equations for and respectively. 

+[x^ + Tfp + Zp - 0 (11) 


l^[(Xpj-Xp)2 + (Ypj-Yy)2] 

p2 

+ ^p - r®-Zp) 


( 12 ) 


Each of equations (11) and (12) admits two solution. Recognizing that the 
cap is totally located above the location Z=0 we choose that X which gives 


"iJ > °- 
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III. 4 Geometric Description of Frequency N Caps 

We are now in a stage where we can qualitatively describe the geometric 
make-up of a cap of any alternate frequency N. Here every pyramid face tri- 
angle has the number of subdivision faces given" by (see, for example 

figure 13 withN=6 for illustration) 

N-1 

= I (2m+l) = (13) 

m=*0 

Hence the total number of faces in the whole cap, S^, is 

Sf = nN^ (14) 

2 

Now since from (13) each of the N subtriangles has three sides and since 
each side is shared by two neighboring triangles, one has the total number, 

Z , of columns to build an isolated single pyramid triangle 

A = ^ (N+1) (15) 

Accordingly, the total number, L, of elements required for the total cap is 
uN 

L = (1+3N) (16) 

Hence, an effective nximber of elements per each of the original pyramid tri- 
angles will be L/n, Pinally, the total number, v, of vertices in the whole 
cap is 

nN 1 

V = [ 1 + (N+1) J (17) 

III. 5- Chord Factors 

Having derived expressions for. the total number of members required to 
build the cap we now indicate that not all of these members are either equal 
or different in lengths. Knowing the minimum number, m, of different lengths 
(also known as the chord factors) is of most importance. Accordingly we now 
proceed to develop formulas f6rsuch a number. Generally speaking, once on 
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FIG. ‘13. NUMBER OF TRIANGULAR SUBDIVISIONS OP 
A PYRAMID FACE (N=6, = 36) 
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the required surface the distances between neighboring vertices will not 
retain their original lengths. This is because of the "stretching" required 
to project them on the surface. In fact, these lengths will be dependent upon 
their original locations on the pyramid face triangle. Various synnnetry relations 
will lead to deriving exact relations for the number of chord factors, m, re- 
quired to build a given antenna. These relations are functions of the frequency 
N and the projection center (Xp, Y^, Z^) but interestingly are independent, 
except for a spherical surface with n=5 (see formulas below) , of n (the number 
of the sides of the original p 3 rramid) . These relations are listed as; 

(a) Vor aparaboloidal surface with projection center (0,0, Zp) where Zp is an 
arbitrary variable 

3 2 1 

m=-^N +N+-^,(N odd, see Fig. 14 for illustrative symmetry) (18) 

3 2 

m‘ = -^ N + N , (N even, see Fig. 15 for illustrative symmetry) . (19) 

(b) For a spherical surface we have the following formulas depending upon the 
projection center and n 

(i) center of projection (0,0,0) and n?^5 

0 O vr Q 

m=-^N odd, see Fig. 16 for illustrative symmetry) (20) 

3 2 N 

m = -^N +’ 2 > even, see Fig. 17 for illustrative symmetry) (21) ‘ 

(ii) center of projection (0,0,0) , n=5 and (Fuller's case) 

1 2 

m = -^ [(N+1) - 1] , (N even and not multiple of 3), Csee Fig, 18 for 

illustrative symmetry) (22) 
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1 2 

m = -^ (N+1) , (H odd and not multiple of 3)5 Xsee Pig. 19 for illustrative 

symmetry) (23) 

m = ^ (^ + 1) , (N multiple of 3) » (see Fig. 20 for illustrative symmetry) 

(24) 

(iii) center of projection (0,0,Zp), with 5 ^ 0 
3 2 1 

m=-^lsi +N+-^, (N odd’, see Fig. 14 for illustrative symmetry) (25) 

3 2 

m-^N +N,(K even, see Fig. 15 for illustrative symmetry) (26) 
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FIG. 14.. NUMBER OP CHORD FACTORS AS PREDICTED 
BY EQUATION (18) 
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FIG. 15. NUMBER OP CHORD FACTORS AS PREDICTED 
BY EQUATION (19) 
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FIG. 16. NUMBER OF CHORD FACTORS AS PREDICTED 
BY EQUATION (20) 
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FIG. 17., NUMBER OF CHORD FACTORS AS PREDICTED 
BY EQUATION (21) 


30 



FIG. 18. NUMBER OF CHORD FACTORS AS PREDICTED BY EQUATION (22) 
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FIG. 19. NUMBER OF CHORD FACTORS AS PREDICTED BY EQUATION (23) 


32 





FIG. 20. NUMBER OF CHORD FACTORS AS PREDICTED BY EQUATION (24) 
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III. 6 Normal Projections from the Intersections of 67 + 60° Array 

An interesting special case of our modeling subdivision and blowing can 
be obtained by choosing thevertexl* whose coordinates are (0,0, R-H) rather 
than the vertex 1 whose coordinates are (0,0,R)for our reference "pyramid". 

If subsequently we project, for any frequency N, from the center (0,0, Z^) 
with Zp -»■ «> we obtain the results for the case where normal projections are 
carried out from each of the intersection of 0, + 60° flat array to meet the 
surface of the solid (see figure 21 for complete -procedure illustration). 

Exact results for this special case can also be obtained as follows: 
by setting = 0, = 0 and = R-H in equation (4) we adapt it to the 0, 

+ 60° flat array subdivision of the aperture circle. For every point (X^jj 
YijjZij) we raise a normal to the flat array which meets the required surface 
at the point which is uniquely determined once the surface is 

specified as follows : For the sphere one has 

and for the paraboloid one has 

Once the points on the surface are specified the calculations follows exactly 
the steps of the alternate breakdown outlined above. 

As will be shown ntimerically later on both the results of the ^act procedure 
of equations (2'3)-and (24) will be' indistinguishable, as far as member lengths 
and smoothness are concerned from the limiting case -of the alternate breakdown 
with the projection center (0,0, Z^) with Z^ It will also be shown that the 

number of chord factors will dramatically decrease for the "normal projection" 
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modeling especially for the case of the paraboloid surface. 

Measure of Smoothness 

As we mentioned earlier the higher the frequency N is, the better the 
discrete system will approximate the required solid surface. An important 
measure of such approximation (smoothness) is the magnitude of the maximum 
distance between the subdivision triangles and solid surface. We refer to this 
measure as 5 . It is obvious that for a given cap geometry the value of <S will 
decrease with increasing frequency. It is also obvious that $ will decrease 
with increasing n. In order to, determine 5 , we proceed as follows'. 

Generally speaking the equation of a flat plane P passing through the 
three arbitrary points and is given by 

= 0 ( 25 ) 

or equivalently by 



AX + BY + CZ + D = 0 
where 



\ "l 1 


^ Zi 1 

A = 

Y2 1 

, B = - 

X2 1 


Y3 Z3 1 


X3 Z3 1 



^1 

^1 

1 


^1 

^1 

"l 

c = 

^2 

^2 

1 

, D = - 

^2 

^2 

^2 



^3 

1 


^3 

^3 

"3 


( 26 a) 


( 26 b, c) 


( 26 d,e) 
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Let us suppose that another curved surface S Cauch as a paraboloid for 
example) passes through these three points. The maximum distance between the 
plane P and the surface S (see figure 22) is the normal distance between the 
plane P. and the tangent plane T to S which is parallel to P; this distance 
is shown as N 2 in figure 22 . The directions ratio of N 2 are given by 

A : B : C (27) 

since N 2 is perpendicular to the planes P and T. Using the above analysis, 
we can determine the maximum derivation for any surface, particularly those of 
the spherical and paraboloidal as follows 
Spherical Surf ace ; 

The equation of spherical surface is given by (see equation 1 ), 

2 2 2 2 

X + Y + Z = R (28) 

Accordingly the direction cosines of any perpendicular to this surface are 
given by: 


8 X ■ 8 Y * 3Z 


(29) 


where the function F is given by 
2 2 2 2 

F = (X + Y + Z - R ) (30) 

At the particular point N 2 (the parallel tangent point), these direction ratios 
must be the same as those of (27) 5 this implies 


8 X * 9Y ’ 9Z 


C 


which from (30) yields 

A : B ; C =X : Y : Z 

Equation (32) can thus be rewritten as 

A = 1 = £ _ 

X Y ^ 


(31) 


(32) 


(33) 
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where X is to be determined. Substituting from (33) into (28) we determine 
the coordinates of the point = (X^s Z^)' as 

X® = A/X , Y® = B/X , 2® = C/X (34) 


and hence X is given by 


X = 


(A^ + + 


ch 


1/2 


/R 


(35) 


Now, the maximum distance 5 between the plane P and the surface 
of the sphere is given by the distance between the point N 2 and the plane P 
which is known as 


g ^ ^fi^:}^Y®iCz5+D 
/A + B + C 
•PaYaboloidal ‘Surface : " 

For the parabolidal surface the function F is given by 


(36) 


2 2 

2 2 ^ 

F = (X + Y + — Z ^)- 

H H 


Using the same analysis used for the sphere one gets 

B ; C 

or equivalently 


3F 

3X 


A 

W + 


(37) 


(38) 


_A = 1_ 

2X 2Y 



(39) 


Accordingly the coordinates of the point are given by 

2 2 

^ - A/C2X) , T® - B/(2X) , z* . E - H((X°) -kA.1 (40) 

where \ is given by 

X = Hc/R^ (41) 

The maximum derivation can again be calculated using the formula 36. 
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•IV; COMPUTER CODE CALCULATIONS 


A general purpose computer code has been written in order to model discrete 

rod-like shells of revolutions as outlined above. This code contains many possible 

combinations of design parameters which we have mentioned in our analysis . The 

code also contains a general and self sufficient graphical subroutine. To see how 

both the computational (modeling) and graphical codes work we have first to define 

their parameters and then proceed to study their flow chart. The graphical subroutine 
is described in Appendix A, 

Input Parameters ; 

F(J>D = focus to aperture diameter ratio . 

AP=D = aperture diameter. 

NS=n = number of equal subdivisions of aperture circle . 

f 1 Sphere 
LT = Shell type: / 

1 2 Paraboloid 

(XP, yP, ZP) = Coordinates of the center of projection. 

/ 1 Single center of projection method 
LN = Design control card: / 

1 2 Normal projection method 

N = frequency of subdivision 

AAl, AA2, AA3 = a^, “ Snler's angles 

EPSL = constant used in the parametric study of chord factors . 

Calculated Parameters : 

R1 = AP/ 2 = radius of the aperture circle 
H = height of the cap 

(XI, Yl, Zl) , (X2, Y2, Z2) , (X3, Y3, Z3) are the coordindates of the represent- 

ative triangle face (1,2,3) (see Sketch 1). 
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1 (XI, Yl, Zl) 



(X2, Y2, Z2) 2 EL2 3 (X3, Y3, Z3) 

Triangular Face 
SKETCH I 

ELI, EL2, a properties of the triangular face. 

X^(r,J), Y^(I,J), Z®(I,J) coordinates of the projected point on the 

surface of the cap. 

Dl(I,J), D2(I,J), D3'(I,J) the side lengths of the sub triangles . 

G(I)=m *= different chord lengths. 

GG(I) = different normalized chord factors. 

GGMX = length of the longest chord, 

GGMN = lenth of the minimum chord. 

IK^ = numerical number of the chord factors. 

KPP=m analytical number of the chord factors. 

LMN = ninnerical number of the parameterized chord factors for a given EPSL. 
GN(LMN) = parameterized chord factors. 

XMAX=6 =» maximum deviation. 
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COMPUTER PROGRAM FLOW CHART 
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WE INTEND TO SUBMIT THE COMPUTER PROGRAM TO 
COSMIC. IN THE MEANWHILE THE PROGRAM WILL 


BE AVAILABLE FROM THE AUTHORS 



V. ILLUSTRATIVE RESULTS 


The utility of our general computer program will be demonstrated by 
generating necessary design parameters of discrete large space antenna re- 
flectors, Specifically we determine the shape and size of flat segmented sur- 
faces which approximate spherical and paraboloidal reflector surfaces. Our 
results will also be displayed in the form of compressions between the spherical 
and paraboloidal designs. Results will be presented for a wide range of pro- 
portions; focus-to-diameter (F/D) radius from .25 to 1 which go all the way 
from a hemisphere to a very shallow dish. Other important parameters which 
may vary are the nximber of pyramid faces, n, the frequency of subdivision, N 
and the aperture diameter D. 

We have carried our numerical calculation on an antenna with the common 

fixed properties D = 100m, n = 6 and N = 10 and the variable properties FOD 

and the center of projection. In Figure 23/ a plot of L and L . for both 

^ max mxn 

spherical and paraboloidal dishes shown as functions of FOB for two differ- - 
ent projection centers namely the center of the sphere (0,0,0) and (0,0, -10 R) 
(this second projection center is equivalent to the normal projection method of 

L^in 

Section 111* 6)..sln Figure 24 various plots of are shown as functions of FOD 

max , 

for the three projection centers (0,0,R— H), (0,0,0) and (0,0, -10 R) . In Figure 

\in 

25 a variation of as a function of projection center is depicted for FOD = 1 

max 

for the paraboloidal surface. As may be seen, the best projection center is that 

L . 

which is around (0,0, -R) where = .985. The corresponding results for 

max 

FOD = 0.25 and for a spherical surface is shown in Figure 26. The variation of 

6 as a function of frequency N is shown in Figure 30 for various FOD values 

for both spherical and paraboloidal surfaces. 

The variation of 6 with FOD is shown in Figure 29 for various projection 

L . 

centers. Finally a variation of as a function of n (for an FOD = 1 

max 

paraboloidal dish) is shown in Figure 28 for two different projection centers. 
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L^in 

As may be easily seen from this figure the can be obtained for n = 6 

Li 

max 

with the projection center being (0,0,0). The variation of this max value 

L . 

of with the frequency is shown in figure 27. Varying the frequency 

max 

will alter very slightly this value. 

From the above figures one can easily draw the conclusion that the 
L 

7 ' ^— ' increases with increasing FOD where both the spherical and the paraboloidal 
mxn 

dishes are practically indistinguishable. On the other hand for lower FOD it 
L 

appears that is much higher for the paraboloid as compared with the sphere. 

min 

Finally, for the plotting illustrations wedepi'cfc in figures 31 and 32 plots 
of a variety of dishes. On each plot we list the necessary parameters used. 
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JJ * 

FIG. 25. VARIATION OF FOR A PARABOLOIDAL DISH AS A 

max 

FUNCTION OF PROJECTION CENTER (FOD = 1.0, n=6, N=l-0 , 
D=100m) 


FIG. 



26. VARIATION OF FOR A SPHERICAL DISH AS A FUNCTION 

■^ax 

OF PROJECTION CENTER (FOD=0.25, n=6, N=10, D=100m) 
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CENTER OF PROJECTION 

( 0 , 0 , 0 ) 

CENTER OF PROJECTION 

y (0,0, R-h) 








0 0^25 0,50 0,75 1^00 


FIG. 29. VARIATION OF NORMALIZED DEVIATION, 6, WITH FOD 
FOR A VARIETY OF PROJECTION CENTERS (n=6 , N=10 
D=100m) 









^^fSaaBwaa^g*^ 



{d) FOD =0.75 



(e) FOD =1.0 



(f) FOD =2.0 


FIG. (31a-f) PLOTS OF PARABOLOIDAL DISHES FOR A VARIETY OF FOD 
(n=6, N=10> D=100m, CENTER OF PROJECTION (0,0,0), 
(a^,a^,a3) = (60°, 30°, 30°) 
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FIG. 32a. VERTICAL VIEW OF A SPHERICAL DISH (n=b, N=50, 
D=725m, FOD = .793 AND CENTER OF PROJECTION 
(0,0,0)) 
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FIG, 32b. SAME AS THAT OP' FIG. 3 2a WITH CENTER 
OF PROJECTION (0,0, R-H) 
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APPENDIX A 


A computer program is written so that we can graphically illustrate 
the cap’s discrete surfaces for any three-dimensional orientation. Due to 
the rotational symmetry about the z-axis we draw all of the cap once we know 
how to draw one of the pyramid's face triangles . 

Let (X^^\ be the coordinates of a point in the first tri- 

angle face of figure Al. The coordinates of the corresponding points on 

( 2 ) ( 2 ) ( 2 ) 

the adjacent faces (counterclockwise) will be , y , z )> 

respectively, (see fig. Al) . 

Notice also 'that Z^^^ = Z • Z'^^^ . Let us use polar coordinates; 

the the point (X^”^\ Z^’^^) can be written in polar form (while suppressing 

the z-coordinate) as 


^(m) ^ ^^(m) ^ ^^i9^i(m-l)(J. 


where m = 1,2,3, ...n 


, 2tt 

(j, = — 


(Al) 


(A2) 


Now suppose that we introduce a new system of axesX, Y and Z (see fig. A2) which 
is fixed in the body of the cap but moves with it.- Also suppose we give to the 
cap any orientation with respect to the fixed system of axes X, Y and Z. 

\ife- must choose the origin, o, of the system X, Y and Z to be the same as the 
origin, ,o of the system Y and z By projecting the cap in the (x-Y) 
plane or the (Y-Z) plane, we obtain a view of the cap in a required orientation. 

The question is knowing the values of X^, Y^ and z^ for a specific point what 
would be the correspgnding values of X^, Y^ and Z^, This can be done by referring 
to the Eulers angles. First let us fix our system of axes in space x»Y and z* 

Then we start our orientation process by letting (x, Y and z) which are fixei 
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in the body of the cap coincide with the fixed system X, Y and Z,. At 

this stage, we start to rotate the body around the Z axes with an angle 

(see figure (A3)). To make it easier to illustrate let us use the notation 

n^ which represent unit vectors along the direction of the fixed axes in the 

-( 1 ) 

body before rotation and let n^ represent unit vectors along the direction 
of the axes after this first rotation. The relation between these unit vectors 
are given by 


— 






cosa^ 

-sxna^ 

n2 


sina^ 

cosa^ 

“3 


0 

0 






^2 

-CD 


(A3) 


r(l) 


Now, let us rotate the body at the new position of X, i.e. around n^ , with 

• -C2) 

an'-- angle to obtain a new set of axes n^ rela^ted to the preceding 

by 


(A4) 

To end our process of orientation, let us rotate the body about the 

-( 2 ) 

new position of Z, i.e. around n^ , with an angle to obtain the final set 
of axes n; \ One notices that this is the set of axes X , Y and Z which 
is fixed in the cap and one also notices that the coordinates of any point 
on the cap ^s surface with respect to -that set of axes are known and related 
to the preceding by 


“1 


^00 


«(2) 

2 

= 

0 COSU2 -sina2 


-(2) 

^2 

•^3 


0 

sxna„ cosa„ 

^ ^ L 


3 
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So, accordingly we finally obtain the relation 


(A5) 





(A6) 


(A7) 


= (sina^CQsa^ + cosa^cosa2Sina^)X^ 

+ (- sina^^ Sina^ + cosa^cosa2COsa2)Y^ 

+ (-cosa^ sina„)Z. (A8) 

^ X Z X 

= (sina2Sina2)S^ + ( sina 200502) + coso2Z^. (A 9 ) 

. So, for different values for the angles a^, 02 and a^, one can obtain different 
view from different angles to the cap. Some results of our plottings are 
illustrated in figures 31.and 32’^ for both spherical and paraboloidal caps. 
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